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INTRODUCTION 

The goal of this research is to develop theoretical, computational, and experimental 
techniques for predicting the effects of irregular topography on long range sound propagation 
in the atmosphere. Irregular topography here is understood to imply a ground surface that 
(1) is not idealizable as being perfectly flat or (2) that is not idealizable as having a constant 
specific acoustic impedance. The interest of this study focuses on circumstances where the 
propagation is similar to what might be expected for noise from low-altitude air vehicles 
flying over suburban or rural terrain, such that rays from the source arrive at angles close 
to grazing incidence. 

The objectives of the project, the experimental facility, and the early progress up 
through August 1987 have been described in the five previous semiannual reports [1,2, 3, 4, 5], 
The present report discusses those activities and developments that have resulted during 
the period, August 1987 through January 1988. 


PERSONNEL 

In addition to A. D. Pierce, Yves H. Berthelot, and G. L. Main, (the three co-principal 
investigators on the project), a visiting scholar, Professor Ji-xun Zhou of the Acoustics 
Institute of the Academy of Sciences of China (Beijing), has been working on the project 
since May 1987. Last fall, the graduate student working on the project, James A. Kearns, 
successfully defended his thesis proposal (a requirement of the Graduate School of Georgia 
Tech). During the past reporting period, Dr. Berthelot, Professor Zhou, and James Kearns 
have been mainly concerned with the the experimental phases of the project, while Allan 
Pierce and Dr. Main have been working primarily on the theoretical aspects. 

Allan Pierce and Yves Berthelot visited NASA Langley Research Center on February 
23, 1987, and discussed complementary NASA and Georgia Tech research activities with 
the NASA technical officer, Dr. John Preisser, with William Wiltshire, and with their 
colleagues. These dicussions resulted in an extension of the original grant for a period of 
one more year (February 1988-January 1989). 
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EXPERIMENTAL VALIDATION OF THE MAE THEORY 

As indicated in the previous technical report [5] on NAG-1-566, a reliable method 
established by Embleton, Piercy and Daigle [6] can now be used to measure in situ the 
acoustic impedance of the diffracting surface used in our scale model experiments. A com- 
parison between theory and experiments was presented in the last technical report [5j only 
for data collected over the high impedance plywood ridge. During the past few months, a 
special effort has been made to analyze thoroughly the more interesting and realistic case 
of the ridge covered with a low impedance material (carpet) with a scaling such that it 
simulates a grass covered hill. The comparison between the results predicted by the MAE 
theory and the experimental data has been presented by Prof. Zhou at the 114 th meeting 
of the Acoustical Society of America, held in Miami, last November. The viewgraphs used 
for the presentation and the corresponding text of the presentation are appended to the 
present report (see p. 41.) The general conclusion of the study is that the asymptotic 
results of the MAE theory (sound field on the diffracting surface, or above it, or far behind 
it) are in excellent agreement with experimental data, and they are computationally simple 
to evaluate. 

Figures 1-3 show some typical results for situations where the ridge is covered with a 
carpet of effective flow resistivity a — 1600 cgs-Rayls. As indicated in the previous report 
[5], the acoustic impedance of the carpet is then determined from the semi- empirical single 
parameter model of Delany and Bazley [7]. Figure 1 shows the sound pressure level (relative 
to the free field expected value) as a function of height above the apex of the ridge, at 5, 10, 
and 15 kHz. The experimental data (dots) follow the predictions of geometrical acoustics 
(solid curve). For comparison, the predictions corresponding to the case of a very high 
impedance surface (plywood) are shown on Fig. 1. 

Figure 2 shows the insertion loss of the ridge covered with the same carpet along the 
surface of the ridge, as a function of a arclength £ nondimensionalized by the frequency 
and the radius of curvature of the ridge. The dimensionless parameter q is a measure of the 
reciprocal of the acoustic impedance normalized by the frequency, the radius of curvature 
of the ridge, and the characteristic impedance of air. (See Reference [3]). Experimental 
data compare well with the predictions of the MAE theory and one can also conclude that 
the method used to estimate the acoustic impedance of the diffracting surface is indeed 
very reasonable, even though the thickness of the carpet is only 6 mm (1/4 of an inch). 
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An interesting result which is predicted correctly by the MAE theory is that, unlike for the 
case of a plywood ridge (high acoustic impedance), the insertion loss on a curved surface of 
low impedance (such as carpet) is frequency dependent. 

Figure 3 shows the insertion loss as a function of nondimensionalized height (tf>) for 
situations where the sound field is mesured far behind (2 meters) the ridge. The dots 
are representative of experimental values while the solid line is obtained from the knife- 
edge diffraction formula, which is an asymtotic version of the more general Fock-Van der 
Pol-Bremmer diffraction integral [3]. The theory predicts that the knife-edge diffraction 
approximation is adequate in the so-called penumbra region (i.e., the transition region 
between the bright zone and the shadow zone). It can be seen from Fig. 3 that, in the 
penumbra region (-1 < ^ < +1) the agreement between theory and experiment is excellent. 
Two interesting results not previously reported in the literature on the subject come out of 
the analysis: (1) the insertion loss along the (slightly bent) axis V* = 0 is correctly predicted 
by the theory to be 6 dB; and (2), the insertion loss far behind the ridge is does not depend 
on the acoustic impedance of the diffracting surface. 
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Sound pressure level as a function of height 
for a given frequency ( X ■ 0 ) 




Figure 1. Sound pressure level as a function of height above the apex of the ridge, at three 
different frequencies (5 kHz, 10 kHz, and 15 kHz). Shown are predicted (solid line) and 
measured (dots) SPL-versus-height curves for sound propagating over carpet (effective flow 
resistivity a = 1600). The dashed line is the prediction for the plywood surface without the 
layer of carpet. 
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Figure 2. Insertion loss on the ridge covered with carpet for four different frequencies (5 
kHz, 10 kHz, 15 kHz, and 20 kHz). Shown are the predicted (solid lines) and the measured 
(dots) SPL-versus-distance £ corresponding to a dimensionless arc length from the apex of 
the ridge. From these data, it is infered that the insertion loss of a diffracting surface of 
low impedance is frequency dependent. 
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Figure 3. Insertion loss in the penumbra region as a function of dimensionless height, 2 
meters behind the ridge, for four different frequencies (10 kHz, 15 kHz, 20 kHz, and 25 
kHz). The penumbra region is centered around a dimensionless height of zero. 
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EXPERIMENTS ON THE APPLICABILITY OF MAE THEORY 


The following pages reprint a paper written and published by Yves H. Berthelot, Allan 
D. Pierce, and J. A. Kearns, during reporting period. The proper citation for this paper is 
as follows: 


Yves H. Berthelot, Allan D. Pierce, and James A. Kearns, “Experiments on the Applicability 
of MAE Techniques for Predicting Sound Diffracting by Irregular Terrains.” paper AIAA- 
87-2668 in Proceedings of the AIAA 11th Aeroacoustics Conference , October 19-21, 1987. 
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Abstract 

The sound field diffracted by a single smooth hill of finite 
impedance is studied both analytically, within the context 
of the theory of Matched Asymptotic Expansions (MAE), 
and experimentally, under laboratory scale modeling con- 
ditions. Special attention is given to the sound field on the 
diffracting surface and throughout the transition region be- 
tween the illuminated and the shadow zones. The MAE 
theory yields integral equations that are amenable to nu- 
merical computations. Experimental results are obtained 
with a spark source producing a pulse of 42 /is duration 
and about 130 Pa at 1 m. The insertion loss of the hill is 
inferred from measurements of the acoustic signals at two 
locations in the field, with subsequent Fourier analysis on 
an IBM PC/AT. In general, experimental results support 
the predictions of the MAE theory, and provide a basis for 
the analysis of more complicated geometries. 

Introduction 

At first glance, it may seem that the propagation of 
sound outdoors is quite simple to analyze. However, this 
superficial attitude changes rapidly when one starts think- 
ing about the challenges arising in real life situations. 1 
For instance, the irregular topography of the terrain may 
cause scattering and diffraction of an incident sound wave. 
The nature of the terrain, i.e., the acoustic impedance 
of the boundary between the terrain and the surround- 
ing medium, may vary with location. At grazing angles 
of incidence, the sound wave may see a non-local reaction 
of the terrain. All the above mentioned difficulties occur 
in practice, even if one neglects the additional complica- 
tions due to atmospheric conditions (temperature gradi- 
ents, wind gradients, turbulence, etc...). The method of 
Matched Asymptotic Expansions 3 - 5 (MAE) is nevertheless 

Copyright © American Institute of Aeronautics and 
Astronautics, Inc., 1917. AH rights reserved. 


a powerful mathematical tool in the analysis of this class 
of problems, provided that the topography varies slowly in 
terms of an acoustic wavelength. In this paper, the diffrac- 
tion of sound by a single curved ridge of finite impedance 
is analyzed with the hope that it may eventually serve as 
a building block for the analysis of more complicated to- 
pographies. Another reason for focusing on the diffraction 
of sound by a single ridge is that there is a direct analogy 
between sound propagation over a flat surface in a curved 
atmosphere (curvature due to temperature and wind gra- 
dients) and sound propagation over a curved surface in 
a fiat atmosphere (with constant properties). Since it is 
much easier to control surface curvature than temperature 
and wind gradients, there is a lot to be gained in modeling 
one by the other. It is therefore crucial to have a good un- 
derstanding of the diffraction of sound by a curved ridge 
under controlled laboratory conditions. The first part of 
the paper deals with the theory of plane wave diffraction 
by curved surfaces of finite impedance, and the following 
section gives a description of the experimental setup used 
to obtain the results shown in the last section. These re- 
sults include essentially the measurement of the insertion 
loss of the ridge as a function of frequency and local radius 
of curvature, both on the diffracting surface and through- 
out the transition region between the illuminated and the 
shadow zones, i.e., in the so-called penumbra region 5 . 

Theory 

The theory of wave diffraction by curved surfaces can 
be traced back to the early 1900’s when scientists were 
studying the transmission of radio waves beyond the line 
of sight and into the electromagnetic shadow of the earth. 
In the 1940’s and 1950’s, V. A. Fock published a number of 
important papers on the theory of the diffraction of elec- 
tromagnetic waves. 5 The following discussion, which has 
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been presented earlier, 6 is largely a translation of Fork’s 
work to the equivalent acoustical problem under investi- 
gation. The derivation, however, is based on the method 
of Matched Asymptotic Expansions (MAE) and it yields 
a uniformly valid solution of the problem. 

Ou ter s olution 

Let a plane wave of constant frequency reflect from 
a curved surface of finite impedance Z,, whose radius of 
curvature R is not necessarily constant, but remains nev- 
ertheless always very large compared to the acoustic wave- 
length. (See Fig. 1). Reasonably far from the surface, in 



where 


\Q- Ui‘/ J rQ- U- f^i . 

F=, + [V1 li T- ' 


with 


and 


Q = (J * 1 + 
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(26) 


(2c) 

(2d) 


The next step is to introduce the following scaled pa- 
rameters L x = R/(kR ) 1 ^ , L v = R/(kR) 2 f* t and Z chMt = 
pc(A:7?) l/3 , so that the functional dependence of the field 
pressure is reduced from five dimensional variables to a 
combination of only three nondimensional variables. In 
other words, one writes the acoustic pressure in the form: 


P 


« P,e’ k 



y_ z \ 

Ly Z chftr / 


( 3 ) 


Since a smooth transition between the solutions far from 
the surface (outer region) and close to the surface (Inner 
region) is required, Eq. (3) is expected to provide an ade- 
quate outer boundary condition of the inner solution. 


Figure 1. Geometry of the problem 


the so-called outer region, the details of the actual wave 
diffraction are overshadowed by the superposition of the 
direct and the reflected waves, so that the field is accu- 
rately predicted by geometrical ray acoustics. Simple geo- 
metrical considerations lead to an expression for the field 
pressure in the form: 

p = p yt ik * + p[A(0)/A(/)| 1/2 (1) 

where P,e'** represents the direct wave, and 9? is the re- 
flection coefficient at the point x 0 ,y 0 on the surface. A(l) 
is a measure of the ray tube area after the wave travels 
a distance l from the reflection point. The resulting field 
pressure can be expressed cxpticitely in terms of the field 
point coordinates (x,y), the acoustic wavenumber fc, the 
radius of curvature R at the apex of the surface, and the 
surface impedance Z # , in the following manner: 

pa P i e ik ’F(x,y,Z„R,k), (2a) 


I nner solution 

In order to get a good picture of the acoustic field 
close to the surface, one has to solve a wave diffraction 
problem, instead of a simplified geometrical acoustics re- 
flection problem. To develop the inner solution, the top of 
the surface is approximated by a parabola y — —z 2 /2R, 
and the Helmholtz equation is expressed in parabolic cylin- 
der coordinates u and v, such that 



X = u[l + («/*)] 

(4a) 


y = u|l + („/2fl)| - u*/2R 

(46) 

so v — 0 
impedance 

corresponds to the diffracting surface, 
boundary condition is also expressed as 

The 


dp 

— + Qp = 0 on the surface, 
on 

( 5 ) 

where Q — tkpc/Z § . The next step is to change from the 
cartesian coordinates x,y to the parabolic cylinder co or- 


dinates £ = u/(2 1 / , t 1 ) and r, = 2 ,/J t >/L„. It can be 
shown 6 ' 7 that the new boundary value problem is now: 
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da <pg 

' "* d 7 t i 


+ = 0 , 


(6) 


with 

dG 

— + qG = 0, at t] = 0, (7) 

where q = i(*/?/2) l/3 (pc/Z,) is a nondimcnsional surface 
admittance. The unknown function £(£, 17 , 9 ) is related to 
the field pressure p by 


P = 


V‘*'»G({,.m). 


( 8 ) 


The outer boundary condition of the inner problem is that 
Eq. (8) matches Eq. (2) at large positive y or large nega- 
tive x. 

The general solution of the above posed boundary 
value problem can be developed by Fourier transform and 
complex variable techniques, with the result 


J [v(a - T)) 

u>; (a) - 9 u>, (a) J ; 

where v(^) and ^(f) are the Fock functions defined by: 


»(f) = w ,/J Ai(f) (10a) 

« v t (f) = 2jr l 'V’' 6 A> (fe^"'' 3 ) (106) 


Eq. (9) is simply related to what Logan 7 terms Fock’s 
form of the van dcr Pol-Brtmmtr diffraction formula. 


can use the wronskian identity w\v - xihv’ = 1 to simplify 
Eq. (9) as follows: 


G(C,0 ,q) 


1 r + o ° e iui 

!-*> «a',(a) - gu),(a) 


( 12 ) 


Note that for 9 = 0, i.e. for the case of a rigid surface, Eq. 
(12) reduces to the special function tabulated by Fock in 
his original paper.* 


The creeping wave series is not absolutely convergent 
at the boundary between the bright zone and the shadow 
zone. In this so-called penumfcra region, it is more appro- 
priate to use an alternative form of Eq. (9). As shown 

by Pierce tt of., 7 the function G, in the penumbra region, 
tends asymptotically, for rj somewhat larger than unity, 
towards a knife-edge diffraction pattern: 

V e * 

- e 1 *’/ 3 !*’'’ + ^—A d (X)](13) 

where X = (2/*)‘/V /4 (£ ~ I 1 ' 3 ) a nd is the 

diffraction integral, which is simply related to Fresnel in- 
tegrals, and which is characteristic of the diffraction by 
sharp edges. The above analysis yields analytical expres- 
sions which are valid throughout the field and which are 
amenable to numerical computations. The following sec- 
tions are devoted to experimental results obtained under 
laboratory conditions, in order to test the applicability 
of MAE techniques for predicting sound diffraction over 
curved surfaces. 


Limiting cases 

Deep into the shadow zone, the diffraction integral 
[Eq. (9)) can be evaluated by a contour deformation and it 
becomes 2t7r times the sum of those residues corresponding 
to poles in the first quadrant, each such term giving rise 
to a creeping wave. The residues series representation in 
the deep shadow zone is given in detail in Ref. 5. 

Another limiting case of interest is the acoustic pres- 
sure on the surface of the ridge, at t] = 0. It can be 
rewritten in the form 

p—P, l e' i ‘ G((,0,q), (11) 

where s denotes the arc distance along the surface from the 
apex of the ridge down the shaded side. In this case, one 


Experimental se tu n 

Laboratory scale experiments are being conducted on 
a 4.88 m by 2.44 m (16 ft by 8 ft) bench made of 2 cm thick 
CDX plywood. The diffracting surface is an arc of cylin- 
der with a radius of curvature of about 2.5 m. The sound 
source is produced by an electric spark across two elec- 
trodes separated by a 1 mm gap, and excited by a voltage 
of 2 kV. The acoustic signature is a repeatable transient 
of 42 ± 2/iS duration, with an amplitude of about 130 Pa 
at 1 m. There are two major reasons for using a transient 
broadband source instead of a continuous single frequency 
signal. First, undesirable echoes from the non-anechoic 
room can be eliminated by properly gating the signals of 
interest; and, second, the frequency dependence of the in- 
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sertion loss of the ridge can be obtained directly over the 
whole frequency range by means of Fourier transform tech- 
niques in a single measurement. Although the spark source 
is much more repeatable than anticipated, the experimen- 
tal procedure does not rely on the consistency between 
consecutive sparks. Instead, two B&K 4136 1/4” micro- 
phones capture the signal at a reference position in the 
free field (unaffected by the ridge) and at a variable field 
point. The ratio of the Fourier transforms of the field and 
reference pressures is independent of the spark signature 
to within 0.5 dB between 8 and 40 kHz. Consequently, the 
ratio of the Fourier transforms, corrected for the difference 
in spherical spreading, the difference in absorption over the 
distances travelled by the signals, and the difference of sen- 
sitivities of the two microphones capturing those signals, 
can be interpreted as a measure of the insertion loss as a 
function of frequency, at a given field point. 

The data acquisition system comprises the two micro- 
phones and their preamplifiers, an analog-to-digital con- 
verter (RC-electronics ISC-16), and an IBM PC/XT with 
a 64 kilobyte memory buffer. All the Fourier analysis 
is performed on the computer by a standard FFT algo- 
rithm, with no window applied to the signals, because it 
was found that spectral leakage was not significant. The 
signals are sampled at a 250 kHz sampling rate, digitized 
with 12 bits resolution, and digitally filtered with a cut-off 
frequency of 75 kHz. 



and comparison with experimental data (•) 

are 2 cm and 0.3 cm respectively, and the horizontal sep- 
aration between the source and the receiver is 100 cm. It 
can be inferred from Fig. 2 that the diffracting surface 
used in our experiment has an effective flow resistivity of 
about 0.8 X 10 6 rayls-cgs. This value is sufficiently high to 
justify the approximation made further in this study, that 
the diffracting surface is perfectly rigid, i.e., that Z, — ► oo 
and g — » 0. 

The first test of the MAE theory presented above is 
shown in Fig. 3 which depicts the insertion loss measured 
on the surface of the ridge as a function of the nondimen- 


Results 


The experimental setup described above is used to an- 
alyze the sound field diffracted by the ridge, both on the 
surface and throughout the transition between the bright 
zone and the shadow zone. For the results reported below, 
the diffracting surface is a 2 cm thick plywood with very 
high impedance Z , . The model of Delany and Bazley* is 
used to determine the impedance of the plywood from the 
value of the effective flow resistance (<r), which is deter- 
mined from a best fit technique described by Embleton et 
al .• Figure 2 shows a set of predicted sound spectra for var- 
ious values of the effective flow resistivity of the plywood, 
together with some experimental data (circles). The ref- 
erence sound pressure level is that which would exist at 
the receiver in the free space. Source and receiver heights 



DIMENSIONLESS ARC LENGTH 


Figure 3. Insertion loss on the surface of the ridge. 
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sional parameter { = {s/R){kRj 2) 1 ^ 3 , where s is the arc 
length measured from the top of the apex of the ridge, R is 
the radius of curvature of the ridge, and k is the wavenum- 
ber. The solid line denotes the predicted dependence (see 
Eq. 12) for the case of a rigid surface; the symbols repre- 
sent data measured at 5, 15, and 25 kHz. The data is in 
good agreement with the theory for positive values of 
(i.e., in the shadow zone). For negative values of £, how- 
ever, the data deviates from the prediction, especially at 
higher frequencies. A possible explanation of this discrep- 
ancy is that the receiving microphone was not perfectly 
flushed with the surface. The prediction obtained from 
the MAE theory is nevertheless quite adequate and it in- 
dicates that the scaling used in the nondimensionalization 
of the problem is correct. 

A second test of the MAE theory is shown in Fig. 
4 which depicts the insertion loss measured at the apex 
of the ridge, as a function of the nondimensional height 
rp = (2kR)(3y/3R) z/7 1 where y is the dimensional height 
above the ridge. The solid curve represents the predicted 
insertion loss (derived from Eq. 2), and the squares and 
circles denote the experimental data at 5 and 15 kHz, re- 
spectively. Both theory and experiment show clearly the 
oscillation in the insertion loss due to the interference be- 
tween the direct wave and the reflected wave. Note that, 
although Eq. 2 is stricly valid only in the outer region 
(large ip) t the agreement between theory and experiment 
is quite acceptable very close to the ridge, as ip — ► 0. 


A third test of the MAE theory is to measure the 
insertion loss in the transition region between the bright 
zone and the deep shadow zone, (i.e., in the penumbra re- 
gion), and compare the results with the knife-edge diffrac- 
tion pattern predicted by Eq. 13. This is the purpose 
of Fig. 5, which shows the insertion loss as a function 
of nondimensional height \p = y{k/nx) l f 2 1 where y is the 
dimensional height above the plywood, and x = 2.09 m is 
the horizontal distance between the apex of the ridge and 
the measuring microphone. The solid line is the result 
of a computation based on Eq. 13, and the squares, the 
empty circles, and the black circles represent data points 
at 5, 15, and 25 kHz, respectively. Again, there is fairly 




Figure 4. Insertion loss above the apex of the ridge 
at 5 kHz (D) and 15 kHz (•) and comparison with 
the theory for q ~ 0. 


Figure 5. Insertion loss 2.09 m behind the ridge. 


good agreement between theory and experiment, and the 
nondimensionalization appears to be adequate. The result 
of a similar experiment is shown in Fig. 6, where the only 
difference with Fig. 5 is the value of x, the distance be- 
hind the apex of the ridge. In Fig. 6, x — 0.64 m, and, 
consequently, the approximation of large ij in Eq. 13 is 
not as good as for Fig. 5, where x = 2.09 m. Also, for 
fairly small values of x, the impedance of the diffracting 
surface is expected to have an effect on the insertion loss, 
and the assumption of a perfectly reflective surface, for 
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the theory shown in Fig. 6, contributes to the discrepancy 
between the theory and the experimental data. Neverthe- 
less, the proper trend is predicted by the MAE theory and 
the nondimensionalization appears to be satisfactory. 



Figure 6. Insertion loss 0.64 m behind the ridge. 


Summary and Conclusions 

The diffraction of sound by a smooth hill has been 
studied both analytically (within the theory of Matched 
Asymptotic Expansions) and experimentally (under labo- 
ratory scale modeling conditions.) The insertion loss at 
any point in the field is expressed in terms of an inte- 
gral called the Fock-van der Pol-Bremmer integral (Eq. 9) 
which involves Airy function of complex argument. The 
applicability of the MAE theory has been checked against 
experimental results in some limiting cases, (on the diffract- 
ing surface, along a vertical axis at the apex of the hill, 
and behind the hill in the penumbra region.) The experi- 
mental results support the predictions of the MAE theory, 
and provide a basis for the analysis of more complicated 
topographies. 
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VIEWGRAPHS USED IN PRESENTATION OF 
THE AIAA-87-2668 PAPER 


The following pages reproduce the viewgraphs that were used by Yves Berthelot in his 
presentation at the 11th aeroacoustics conference in Sunnyvale, Ca (October 19-21,1987) of 
the paper 

“Experiments on the Applicability of MAE Techniques for Predicting Sound Diffracting 
by Irregular Terrains.” paper AIAA-87-2668 in Proceedings of the AIAA 11th Aeroacoustics 
Conference, October 19-21, 1987. 


The text of this paper appears in the preceding section. 
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OUTLINE 

0. Introduction 

1. Theoretical model 

2. Experimental method 

3. Surface impedance at grazing angles 

4. Experimental results 


5. Conclusions 
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MOTIVATION : PREDICT THE SOUND FIELD OF LOW- 
ALTITUDE FLYING AIRCRAFTS . 





irregular topography 


METHOD : MATCHED ASYMPTOTIC EXPANSIONS (MAE) 

m 

PROTOTYPE PROBLEM: SINGLE CURVED RIDGE 



BRIGHT SON! 
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OUTER SOLUTION: GEOMETRICAL ACOUSTICS 



In the vicinity of the top of the ridge: 


P* Pi e tt * F(x,y y Z,R y k) 


With seating: 

L. = R/(kR)''> ; L, = R/(kR)*l> ; Z^. - fc(kR)'/> 
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INNER SOLUTION: wave dijfractioa 


{ V’p+tfpm 0 

~ + Qp * 0 on the waibct 

with Q « ikpe/Z 

Change to parabolic cylinder coordinates rj 


if *'’- 0 


A 


at *»• 


with: 


P* e< ** e**'* G(t t ij,q) 




P£ 
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GENERAL SOLUTION | 

MATCHING PRINCIPLE : AT LARGE rj THE INNER SOLUTION 
MUST MATCH THE OUTER SOLUTION. 





®(x) * » l/ *Ai(x) 

Wi(x) m t*"** Aifx* 4 *^ 1 ) 


FOCK - VAN DER POL • BREMMER DIFFRACTION INTEGRAL 
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LIMITING CASES 


1) Shadow aona (( > rp!*) 


0{( t if t q) • H* £(i**Mm* a* p*k* el fctfapaod) 

It reduce* to the remit* * Keller (1956) and Hap* et mL (1978). 
9) Penumbra (( ~ rft*) 


G ps *-«*/» e*» - [HiX^-**™** + 








THE RATIO n(itf)/p,(h>) IS INDEPENDENT 
OF THE SPARK SIGNATURE . 

IT IS THEREFORE A MEASURE OF 
THE INSERTION LOSS OF THE RIDGE. 







SIGNAL PROCESSING 
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/ 
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. r 77 


7 


7 / / 


/ ' 


• For two different surfaces 

(a) Hard ( 1/4 ” thick plywood) 

(b) Soft ( 1/4 ” thick carpet) 


• Measure the insertion loss 

(a) vertically, at the apex 

(b) vertically, behind the ridge 

(c) vertically, far behind the ridge 

(d) on the surface of the ridge 
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3. Surface Impedance 



Ref: T.F.W. Embleton, et a I. JASA 74 (4) 1983, 1239-1244. 



a = effective flow resistivity 




SOUND PRESSURE LEVEL - dB 
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(0 



('Pl^ood ) 



SOUND PRESSURE LEVEL 
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FREQUENCY (Ha$ 
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GRAZING ANGLE 
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Scaling 


• Experimental results: 

- Frequency range 4.0 - 20 kHz 

- Grazing angles 1° - 15° 

Surface can be treated as locally reacting 
with an effective flow resistivity <x 

-a — 1,600 cgs-Rayls (for carpet) 

- er = 80,000 cgs-Rayls (for plywood) 

• Scaling: 400 - 2,000 Hz 

- Gear pet = 160 cgs-Rayls — ► grass 1 ** 

- <?pi V wood = 8,000 cgs-Rayls — ► hard packed earth 1 ** 


(,) : From T.F.W. Emblet on, et al., JASA 74(4), 1983, 1239-1244. 



RELATIVE SOUND LEVEL • 
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XNseeTion toss 

A&we THE 





HEKHT(m) 


tWot^ (terpe - ) 



JSERTION LOSS ~ dB 



NAG-1-566 Semiannual Report No. 6, page 37 





DIMENSIONLESS HEIGHT 


INSERTION LOSS 



INSERTION LOSS 





NAG- 1-566 Semiannual Report No. 6, page 40 


CONCLUSION 


Laboratory scale experiments indicate that the MAE 
method is adequate to predict the sound field diffracted 
by a curved boundary of finite impedance. 


The Fock-Van der Pol-Bremmer diffraction integral can 
be efficiently computed at any field point without hav- 
ing to solve for the whole field (e.g., finite elements). 
It also reduces to some simple form for some special 
regions of the field (e.g, outer region, far behind the 
ridge, deep in the shadow zone, and on the diffracting 
surface). 


The MAE method can be used as a “building block” 
for the analysis of more complicated geometries (e.g., 
undulating terrain). 
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THE EFFECT OF FINITE SURFACE ACOUSTIC IMPEDANCE 


The following pages reprint a paper that has been presented at the 114th meeting of 
the Acoustical Society of America, held in Miami, November 16-20,1987. Professor Zhou 
gave the paper. The proper citation for this paper is as follows. 

Ji-xun Zhou, James A. Kearns, Yves H. Berthelot, and Allan D. Pierce, “The effect of 
finite surface acoustic impedance on sound fields near a smooth diffracting ridge,” presented 
at the 114 th meeting of the Acoustical Society of America, Miami, Fa (1987). Abstract in 
J. Acouat. Soc. Am. Suppl. 1, 82, S49-50 (1987). 
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The effect of finite surface acoustic impedance 
on sound field near a smooth diffraction ridge 

Ji-xun Zhou 
James A. Kearns 
Yves H. Berthelot 
Allan D. Pierce 

(School of Mechanical Engineering, Georgia Tech) 


I. Introduction 

II. Summary of the MAE diffraction theory 
by a smooth ridge 

III. Experimental apparatus 

IV. Experimental results 

V. Some conclusion 


Slide 1 
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(I) Introduction 



Sound Propagation over undulating terrain 



Simplified Diffraction Model 


Slide 2 



NAG- 1-566 Semiannual Report No. 6, page 44 


(II) The summary of the MAE diffraction theory 
by a smooth ridge 

(2,1) Outer ray resolution 

( in the vicinity of ridge’s top ^ 



^ = (2Jt//2 3 )[<? 3 - (8/27 )* 3 - (2/3)JRxy] 


Slide 3 
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(2.2) General Inner wave resolution 


p = Pie ik * e * 1 ' 3 G((, V ,q) 


g ^V,9) = -fa £ [«(« - V) ~ 4(a) ) -wl) U, ‘ (a " ’ )) ] e ‘° t<ia 


v(x) = * 1/2 Ai(x) 
un( x ) = e"'* 2 k 1 12 Ai(xe iw/2 ) 

q = i(kR/2) l / 3 pc/Z s 

FOCK - VAN DER POL - BREMMER DIFFRACTION INTEGRAL 


Slide 4 
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(2.3) The sound pressure on the surface of a ridge 


roo 

G{Z,O t q) = 7T“ 1/2 / — 

» o 


2 [0) - e i2n / 3 qw 2 (P) 

• »< 


dp 


roo .•#( 

+ *"‘ /2 / -T7-1 5 

Jo W 1 («) - 9«>i ( Q ) 


w » ( z ) — 2x 1 / 2 A.\{zc i7n / 3 ) u>[ (*) = 

w 2 (r) = e~ in l*2n 117 A\{ze~ i7 "l*) ti4(z) — 



e «5w/e2 Jr l/2AJ , ( ze *'2ir/3) 

e -ia»/« 2jr i/a A i'( ze -i2*/3) 
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(2.4) The sound pressure in the transition region 

P =e -.r/v<- 

m 

where X ^(2/t) 1 / V /4 (f-»7 l/2 ) 

y 

S ^ = (u|R)(kR|t)' l, _ 

x penumbra 

>7 = (»/«) ( *’ R> /2) ' ' 

s 77 V >w 

M = 7 ^ [- (2fly + tf 2 ) + v/ (272® + )» + 4i2 3 ] ‘ ^ 

v = i? [(x/u) - l] 

= 7-Rn. at £ = rj 1/2 


Slide 6 
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(III) Experimental apparatus 



Block diagram of the experiment. Spark source (S) 
includes power supply (P), resistor (R), and capacitor (C): 
signals from reference microphone (A/j) and field microphone 
{M2) are amplified (A), digitized (A/D), and processed by the 
IBM PC/XT. 


Slide 7 
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(IV) Experimental results 

(4.1) Surface acoustic impedance 

P_ _ ( e ik ' r ' \ ( e ik ^ \ 

Po V k x r x ) p \ k x r 2 ) k x r 2 

F(u>) = 1 + iir Ui w in e ~ w erfc( - iyfui) w ~ i ^ ~ cos 2 <f> J 


Earlier work by Embleton» Piercy and Daigle 


Impedance model of Delany and Bazley 

0.73 



r/ P,c ‘ 

X 2 /p jC, = -f 1 1 .9 (f /a) ~ 0 73 Zj 

a 2 /k t = 1 + 10.8 \f/a)~ 010 
0 2 /k t = 10.3 (/*/ a )- 039 
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SOUND PRESSURE LEVEL - dB SOUND PRESSURE 










RELATIVE SOUND 
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Relative sound level as a function of frequency 
for a given height ( X * 0 ) 



Slide li 
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CMN3CKISS ARC LDCIN 


Insertion loss along the ridge covered with carpet 
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7TT 


•% 






r r r 

ran 

height ( X - 2m ) 
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(V) Some conclusions 

1. The comparison of the MAE diffraction theory 
with experiment from a smooth ridge is excellent. 

2. For small grazing angles and the experimental 
frequency range, the carpet backed with the plywood 
can be treated as a local reaction impedance. 

3. The insertion loss along a finite impedance ridge 
is more sensitive to the frequency than the 
insertion loss along a rigid ridge. 

4. In the penumbra, far behind the ridge, the 
diffraction of a smooth ridge with a finite 
impedance is the same as a rigid knife-edge diffraction. 


Slide 18 
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THESIS PROPOSAL OF JAMES A. KEARNS 


The following pages reprint the thesis proposal successfully defended by James A. 
Kearns on November 12, 1987. The thesis committee was made of Drs. Berthelot (Chair- 
man), Pierce, Main, Jarzynski, Zhou, Long, and Benkeser. (The thesis proposal is a require- 
ment of the Graduate School of the Georgia Institute of Technology which is an intermediate 
step toward the completion of the doctoral degree.) 
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INTRODUCTION 

The principal goal of this research is to gather and analyze experimental data that is useful 
in understanding the effects of irregular topography on long range sound propagation in the 
atmosphere. The experiments are focused on situations where the rays associated with the 
propagating sound wave are impinging upon the topography at angles near grazing incidence. 
Such a situation may be visualized as that when a low flying aircraft produces an acoustical 
disturbance over the nearby terrain. An analogous situation arises when wind or temperature 
gradients result in curved propagation paths in the atmosphere. In either situation, the interaction 
between the impinging sound field and the surface topography is expected to be quite complex. 
In light of this expectation, a much simpler model will be examined. The model chosen is that of 
a single, cylindrical ridge lying on a flat table top. The expectation is that a complex topography 
may be theoretically reconstructed as a set of adjoining single ridges. A spark generator provides 
an acoustical point source. Experiments will be conducted to describe the field on, above, and 
behind the ridge. The data will be compared with numerical predictions derived from the theory 
as presented by Pierce. 1 

This proposal consists of three major parts. The first of these is a synopsis of the rele- 
vant theory and its history. This includes a survey of influential papers since 1945. The second 
part is a review of previous and relevant experimental work unhiding studies of diffraction over 
cones and other bodies of revolution. Finally, the third part contains descriptions of the proposed 
experimental and computational work. The proposed computational work is presented immedi- 
ately after the section on the history of the theory. Likewise, the proposed experimental work is 
presented immediately following the review of previous experimental work. A summary of the 
objectives of this work is presented at the end of the text. 

Dr. Ji-xun Zhou has been working and will continue to work closely with the author. He 
has assisted in all aspects of the project including the capturing of data and the computation of 
theoretical results. 
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SURVEY OF RELATED THEORETICAL WORK 

The study of high frequency diffraction over curved surfaces dates at least as far hack as the 
most recent turn of the century. At that time there was a good deal of interest in studying the way 
in which radio waves propagated over the surface of the earth. A typical model of this phenomenon 
was that of a vertical electric dipole in the presence of a metallic sphere. In all cases, the radius of 
the sphere was assumed to be much greater than the wavelength of propagation. Unfortunately, 
under such circumstances the eigenfunction expansion of the solution in the shadow’ zone is slowly 
convergent. However, Watson (1918), w*ho drew upon the work of Poincare (1910) and Nicholson 
(1910), managed to express the expansion in terms of a residue series. From this residue series, 
an associated contour integral was constructed. Watson showed by deforming the contour in a 
particular way that a new, faster converging resudue series can be derived. This process is known 
as Watson’s transformation and is still used to evaluate the solution deep in the shadow zone. 
At points near the geometrical shadow boundary, this series solution is often slow to converge. 
WTiite (1922) extended the work of W 7 atson for cases when the observer is outside the shadow 
zone. 2 

In the late 1930's, van der Pol and Bremmer (1936-38) published several notable papers on 
"he propagation of radio waves over the surface of the earth. Shortly thereafter, Fock (1945) 
published a paper in which he put forth the concept of the local field in the penumbra region. 
Fock considered the case of a high frequency wave incident upon a perfectly conducting surface 
of continuously varying curvature. By way of a physical argument, Fock claimed that the surface 
current distribution in and near the shadow' boundary depends upon the local curvature and upon 
the incident field. This dependence was expressed in the form of a “universal” function 

G(0 = [ -^- da (1) 

Jr <( Q ) 

where 1) f is a representative distance of a given point on the surface to the geometrical shadow* 
boundary, and 2) T is the contour which runs from infinity to the origin along the path Arg z — 
2r /3 and from the origin to infinity along the real axis, (a) is an Airy function and is a bounded 
solution to tt* M (cr) = ow(a). The function G(£) decreases exponentially with f for large positive 
values of f in the shadow zone and asymptotically approaches 2 for large negative values of £ in 
the illuminated region. The characteristic width of the penumbra was given as d = [^o 2 ] 1 / 3 * 
w'here R 0 is the radius of curvature of *he surface shadow* boundary in the plane of incidence. In 
the following year, 1946, Fock derived the same result directly from the reduced wave equation. 
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In addition, Fock extended the solution to the region in the “neighborhood” of the surface and 
included bodies which are “good” conductors in the sense of the Leontovich boundary condition. 3 

Approximately one decade after the publication of Fock’s work, Keller (195C) published the 
“Geometrical Theory of Diffraction.” Keller’s theory was based upon a heuristic application of 
Fermat’s principle. In order to account for diffraction effects, Keller added a modified form of the 
principle which states that “the diffracted rays connecting the points P and Q are those curves 
which render stationary the Fermat integral among all curves connecting P and Q, and having 
a point or arc on the diffracting obstacle.” 4 (Note: The original wording was arc on the 
cylinder.” since this was the geometry originally examined.) This statement can be aplied to any 
body of arbitrary shape. The resulting solution is in terms of unknown “diffraction coefficients”, 
which are analogous to reflection coefficients, and which are found by comparison to well known 
exact solutions of various ranonical problems. These coefficients are functions of the incident 
field, the wavelength, and the surface properties and geometry at the point of diffraction. These 
coefficients follow the spirit of Fock’s notion of the local field in the penumbra. Keller’s theory 
formalized the concept of “creeping” waves, so-called by Franz and Depperman (1952), for an 
arbitrarily shaped body. The “creeping” wave hypothesis is that a surface wave is produced 
at the point of diffraction and that this wave proceeds to travel into the shadow zone along a 
geodesic of the surface. This wave continuously sheds rays tangential to every point along the 
geodesic in accord with the modified Fermat’s principle of diffraction. Keller’s theory is a high 
frequency approximation although good results have been reported at wavelengths on the order 
of the dimensions of the diffracting obstacle. This theory is not valid near the geometric shadow' 
boundary or on the surface of the diffracting obstacle. 

In 1959, Goodrich examined Fock’s theory in light of the work of Keller, Franz and Dep- 
perman, and Kazarinoff and Ritt (1959). Goodrich showed that by changing Fock’s choice of 
coordinate transformation that Fock’s theory .will yield properly the creeping wave behavior dis- 
cussed above. Also, Goodrich noted that Fock’s theory is essentially a 2-dimensional one in that 
it does not account for ray convergence on the surface of a 3-dimensional body. This difficulty 
was overcome in the derivation of Hong (1967). 6 

A respected analysis of previously published diffraction theories was made by Logan (1959). 
Subsequently, Logan and Yee (1962) presented a thorough mathematical treatment of the inter- 
relationships between the relevant theories for diffraction from a convex body. More recently, 
Ivanov (1971) examined plane wave diffraction from an ideally reflecting cylinder. He proceeded 
to match an asymptotic solution valid in the shadow' zone to an asymptotic solution valid in the 
illuminated region. Ivanov’s matched solution w f a$ continuous but not necessarily smooth at the 
shadow' boundary. Later, Pathak (1979) found his own uniform, asymptotic solution. His solu- 
tion is for a perfectly conducting surface and is valid everywhere e ~ept at the surface. Pathak’s 
solution is given in terms of well tabulated functions similar to those of Fock. 
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The current theory of interest is that as presented by Pierce (1980). His solution employs 
the method of Matched Asymptotic Expansions to derive a uniform solution that 1) is valid on 
much of the surface, and 2) reduces asymptotically to the geometrical acoustics solution in the 
illuminated region. A diagram of the situation and the associated regions of interest is shown in 
Fig. (1). The general solution has the form 

p = P t -e tt V«*/» <?(*,,,,) (2) 

and 

= '" /! /_ !" (o - *> - ' ,, ] £ ‘“ da 131 

where r(f) apd itfj(C) are “Fock functions” and simply related to Airy functions of complex 
argument* £ and q are nondimensionalized and transformed coordinates, and g is a nondimen- 
sionalized surface admittance. The variable it is proportional to £. The above integral solution 
(3) is termed by Pierce the Fock-van-dcr-PoI-Bremmer function and is “trivially related to what 
Logan (1959) terms Fock‘s form of the van der PoUBremmer diffraction formula”* 

In the limit of an acoustically ridge surface, which is analogous to a perfectly conducting 
surface for cases of electromagnetic propagation, the above integral function, (?, reduces to 

( 4 ) 

on the surface of the ridge. This result is nearly identical to that derived by Fock (Eq. 1). Fock 
includes the factor e'* ^ as part of the function G while Pierce separates it from G as in Eq. (2). 
The differing contour paths are the only real discrepancy. 

Proposed Computational Work 

The proposed computational work consists of finding solutions to Eq. (3) for various ranges 
of the parameters £ and 17. Algorithms for many of the ranges of interest have been previously 
derived by Pierce and Main (1987). These algorithms will be used and their solutions compared 
with the experimental data. New or improved versions of these algorithms will be applied to those 
cases in which significant discrepancies between theory and experiment are found. A program 
for the evaluation of the Airy function for arbitrary complex argument has already been written 
by Pierce. This program will be utilized as a base for the other computations. Any necessary 
refinements to this program will be carried out. The results will be compared against experimental 
data for certain regions of space on, above, and behind the ridge. The particular regions are 
outlined in the section on proposed experimental work. 



\y Nr Y V Y h/Nk > V 
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1. The different regimes behind and near the top of a curved surface when a plane 


wave is incident. 
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SURVEY OF RELEVANT EXPERIMENTAL WORK 

A survey of related experimental studies is presented in this section. This includes exper- 
imental studies of the diffraction of waves incident upon cones and other bodies of revolution 
as well as experiments on wave propagation over and on curved surfaces. A significant portion 
of the relevant experimental work has been carried out in the context of electro- magnetic wave 
propagation. An effort is made to include this work although an emphasis is placed upon work 
performed in the realm of acoustics. To the best of the author's knowledge, the survey includes 
most of the relevant experimental work from the past two decades. The goal of this survey is to 
summarize the recent experimental work. 

Cones and other bodies of revolution offer simple 3-dimensional geometries which 1) are easy 
to find or construct and, thus, allow for experimental testing of theories such as the Geomet- 
rical Theory of Diffraction (GTD), and 2) readily admit numerical solutions due to symmetry 
considerations. For this reason, many recent papers have reported the use of cones, spheres, 
and ellipsoids in their model experiments. Keller (1961) made a comparison between theoretical 
and experimental results for the backs cat tering cross-section produced by the electromagnetic 
irradiation of a perfectly conducting cone. The experimental results were obtained by Keys and 
PrimHi (1959) for six different cones with half-angles ranging from 4 to 20 degrees. The data is in 
good agreement with predictions based upon Keller’s Geometrical Theory of Diffraction. Bechtel 
(1965) carried out a series of experiments to measure the radar cross-section (RCS) of cones with 
half- angles of 4, 15, and 90 degrees. The size of the base, a, was also varied between 0.98 and 
2.87 wavelengths in diameter. Again, Bechtel compared his results with predictions based upon 
Keller’s GTD, For ka on the order of 8 or 9, Bechtel found good agreement between the theory 
and experiment except for cases when the cone is observed within 30-40 degrees of nose-on with 
a vertically polarized wave (i.e. wavefront parallel to the axis of the cone). For smaller cone 
bases on the order of ka = 3, the predicted RCS matched the observed RCS to within 5 dB but 
the predicted shape or form of the RCS was not well observed, particularly within 30 degrees 
of nose-on. Another set of reported data on cone diffraction was gathered by Bargeliotes et al . 
(1975). Mittra and Safavi-naini (1978) compared these experimental results with theory which 
considered the potential field produced by point sources on the surface of a cone. The array of 
point sources represented the surface currents as derived asymptotically by Fock and others. The 
comparison showed good agreement when the diffracted rays had traveled a significant angular 
distance around the cone. No result s for * ar grazing angles of incidence were found in their 
report. An extensive study of acoustic scattering from various bodies was made by Lang (1980). 
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He measured forward and back scattering as well as creeping wave type phenomena from smooth, 
rigid objects such as baffles, cylinders, cones, and cone variants. 

Experimental data on high frequency diffraction and/or reflection from other bodies of rev- 
olution has been reported by a number of sources. Neubauer (1967) measured the travel time 
from source to receiver of a sonic pulse diffracted around a circular aluminum cylinder. The 
measurements were made in water. He positioned a narrow-beam source such that the central 
ray arrived tangential to the cylinder surface and normal to the axis of the cylinder. An array 
of five equally spaced hydrophones measured the amplitude and arrival time of the pulse. The 
hydrophones were deep in the shadow zone, Neubauer found that his results were in excellent 
agreement with those predicted by the “creeping wave” theory of Franz and Keller. Neubauer 
also reported that a slit of less than a wavelength in width and positioned along the grazing line 
allowed a “...large part of the wave to pass.” Further, he placed a baffle against the cylinder 
surface to block passage through the water and, subsequently, found little or no transmission. 
Thus, he concluded that* the wave was transmitted on the water side of the cylinder boundary. 
Foxwell (1970) and Blake and Wilson (1977) carried out related experiments. Foxwell measured 
the diffracted field on the shadowed surface of a rigid sphere. A point source located on the 
surface formed the illuminated pole. At high wavenumbers, Foxwell observed interference effects 
predicted by creeping wave theory. Blake and Wilson performed an analogous experiment using a 
highly eccentric prolate spheroid. Their sound source was located at a distance from the surface 
equal to the length of the spheroid and along the major axis. Measurements revealed an illumi- 
nated spot at the antipole and the existence of the creeping wave interference pattern near the 
antipole. In addition, Blake and Wilson report that 1) the measured levels at the antipole “de- 
crease roughly as the reciprocal of frequency”, 2) that the shadow zone measurements are “well 
approximated by Keller’s Geometric Theory of Diffraction”, and 3) that deviations on the order 
of 3 dB from the theory “are observed at high frequencies and at coordinates off the antipole.”. 
Lang (1980) examined acoustic scattering from the surface of a thin prolate spheroid and used the 
results to deduce characteristics of acoustic scattering from any smoothly curved surface. Lang 
concluded that the “backscattered pressure from smoothly curved bodies is determined almost 
exclusively by specular effects, even at wavelengths that are relatively large compared to the ap- 
propriate dimensions of the scatterers.” He also observed “evidence of surface fields very close to 
smoothly curved scatterers” that he “identified with Franz-type creeping waves.” Finally, Lang 
derived the Freedman theory of echo formation and compared the predictions of this theory with 
experimental results. He found that the theory “produces a reasonably good model of scattering 
from baffles but . . . exhibits serious errors in cases involving three-dimensional bodies.” 

More recently, Almgren (1986) measured the insertion loss above and behind convex and 
concave cylinders. A distant spark source provided a near planar incident wave. Almgren was 
interested in examining the analogous relationship between sound propagat’^n over a curved 
surface and that in a medium with a vertical sound speed gradient. He found that “it is reasonable 
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to simulate the effect of refraction due to a sound speed gradient ... with a curved ground scale 
model.” Only small errors on the order of 3 dB or less were found for all measurements including 
those at grazing incidence. These errors were relative to the theory as presented in the works 
of Pekeris (1946), Pridmore and Brown (1962), and Pierce (1981). Berry and Daigle (1987) also 
were interested in the refraction of sound in a stratified atmosphere. Like Almgren, they chose to 
examine the analogous case of diffraction over a convex curved surface. A tone burst mechanism 
was used to produce a narrow band point source. Their measurements were confined to listener 
positions deep in the shadow zone and in the penumbra region while the point source was placed 
either on or slightly above the surface of the ridge. Measurements were made over a frequency 
range of 0.3 to 10 kHz. For listener positions below the shadow boundary, they compared their 
results to predictions based upon an extended version of the creeping wave theory. For listener 
positions above the shadow boundary, they compared their results to predictions based upon 
geometrical theory. At the shadow boundary, they reported that the two theories agreed to within 
\ dB. However, the experimental results differed from the theoretical results by as much as 5 dB 
when they were compared for listener positions on or near the shadow boundary. Comparisons 
outside of this region were good. 
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PROPOSED EXPERIMENTAL WORK 

The intent of the proposed work is to investigate the acoustic scattering at large angles 
of incidence from a curved surface, with particular attention paid to transition phenomena. A 
laboratory and experiments have been designed in order to make such an investigation. The 
experiments can be grouped into four classes. They are: 

1) Preliminary experiments to determine “free field” propagation properties such as the 
sound speed, source directionality, attenuation, and non-linear steepening of waveforms; 

2) Experiments to determine the specific acoustic impedance of two or more surfaces; 

3) Experiments to determine the insertion loss on, above, and behind the model curved 
surface; and, 

4) Experiments to determine the effect of echoes from one ridge on the insertion loss near 
another ridge for a series of two ridges aligned back to back. 

The majority of the proposed work consists of carrying out the above experiments and comparing 
the acquired data to the results predicted by the theory presented earlier. 

Experimental Facility 

The following is a brief description of the laboratory and some common components of the 
experimental methods. 

Four tables were made to be used for the scale model experiments. Each table is 1.2 meter 
wide by 2.4 meter long and 0.9 meter high. These four tables were bolted together, forming one 
table 4.9 meter long by 2.4 meter wide. The table top is CDX plywood, 2 cm thick. A curved 
surface to be mounted on the table was constructed of exterior plywood and used as a laboratory 
scale model of a topographical ridge. The contour of the surface has the shape of an arc of a 
circle wdth an approximate radius of curvature of 2.5 m. The complete laboratory including the 
table with the curved surface (ridge) mounted upon it is depicted in Fig. (2). 

A sound source is produced by an electric spark generator. The design of the spark generator 
is a modification of a design by Dr. Mendel Kleiner at the Chalmers Institute of Technology. 
The spark gap is approximately 1 to 2 millimeters wide. A gap voltage of 2 kV is sufficient for 
discharge. Typical peak acoustic pressures at a distance of 1 meter are on the order of 120 to 130 
Pascals. Typical rise time and characteristic pulse duration are of the order of 10 /is and 42 /is, 
respectively. The resulting sound source is broadband r 2 - 40 kHz) with the peak level residing 
in the range of 20 - 25 kHz. 
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The data acquisition system is composed of two Bruel k Kjaer quarter-inch condenser mi- 
crophones, two amplifiers, an analog- to- digital converter, and an IBM personal computer. The 
amplified analog signal of each microphone is converted to digital form by an integrated hard- 
ware and software system produced by RC Electronics Inc. and called “Computerscope ISC-16." 
This system consists of a 16 channel A/D board which is inserted into the IBM PC, an external 
instrument interface, and the scope driver software. The system has a 64 kilobyte memory buffer 
and is capable of gathering data at an aggregate rate of 500 kHz. 

General Procedure 

The general experimental procedure is similar for all of the experiments described heie. Most 
of the experiments use two microphones, one at a reference point in the free field, and the other 
at a desired field point. The two corresponding analog voltage signals, which are output from the 
microphones, are amplified, sampled at intervals of 4.0 ^s, and digitized to 12 bit precision. The 
sampling process is triggered such that the source waveform arrives at each microphone within 
the limited time window afforded by the finite memory buffer. A constant is automatically added 
to each data sample such that the mean of the data sample is approximately zero. Since the 
voltage increment registered by a microphone is opposite in sign to that of the corresponding 
pressure increment, the sign of the digitized data is reversed. Then the data sets are filtered by 
a lowpass digital filter with an upper cutoff frequency of 75 kHz. Waveform portions of interest, 
and representative of what would be received if there were no undesired reflections contaminating 
the data, are selected. Each portion is centered in a field of zeroes, and thus two new abbreviated 
waveform windows are created. These windows have the same duration (usually, 500 times 4.0 
/is), but usually not the same time beginnings. The replacement of the extra data points by 
zeroes is in accord with the expectation that an acoustic pulse has negligible residual effect on 
the ambient pressure, and that the acoustic pressure before and after each waveform should 
ideally be zero. The digital Fourier transform of each window is computed and interpreted as if 
the value of the acoustic pressure were identically zero before the start and after the end of each 
window. 

Preliminary Experiments 

A set of preliminary experiments have been performed in order to establish 1) the ambient 
speed of sound, 2) the directionality of the source, 3) the magnitude of propagation losses due to 
spreading and absorption, and 4) the extent of non-linear effects. The experiments were performed 
over the flat table top, with a fixed source height equal to the height of the ridge, and with a 
common reference microphone position. Results show that the ratio of the Fourier transform of 
the field signal to that of the reference signal is generally independent of the particular source 
signal. Also, the source is within ± 1 dB of being omnidirectional over a range of -12 to 24 de- 
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2. Schematic of the interior of the laboratory room used in the study, showing the 
sequence of data processing associated with a single firing of the spark source. 
Reference and field microphones are shown stationed in front and behind of 
the ridge, respectively. A triggering device activates the data capturing sys- 
tem simultaneous to the firing of the spark source. The analog signals of the 
two microphones are sampled at a rate of 250 kHz, digitized to 12 bit preci- 
sion, stored in a 64 kilobyte memory buffer, and low pass filtered. Interesting 
portions of the total data field are then input to a discrete Fourier transform 
program. The output of this program is interpreted as the Fourier components 
of the original waveform portions. 
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grees from the horizontal. The local speed of sound was deduced after measuring the arrival time 
of a pressure wave at 11 different downstream positions. The measurements were made at equal 
intervals over a downstream range of 1 to 3.5 meters. 

The same set of 11 data points were used to deduce the free field dependence of the ratio of 
the Fourier transform of the pressure waveform at the field microphone to that at the reference 
microphone. The results show that free field propagation losses can be approximated by those 
attributed to spherical spreading and absorption in air. Finally, the shape of the pressure waves 
were observed for signs of wave steepening while the magnitude of their Fourier transforms were 
observed for redistribution of acoustic energy from lower frequencies into higher frequencies. 

Specific Acoustic Impedance Measurement 

Several experimental methods were applied to determine the acoustic impedance of both the 
bare plywood table top and the table top covered with a low-cut commercial carpet. Later, the 
acoustic impedance of the surface of the ridge will be assumed to be roughly the same as that of 
the table top. The goal of these experiments is to develop and/or test a procedure for acquiring 
a reasonable estimate of the acoustic impedance of a surface at grazing angles of incidence. 

The first method is referred to as the “Direct” method. In this method, the pressure pulse 
arrives at a single microphone along two paths, one direct (distance R \ ) and one reflected (distance 
Rz ). The paths are chosen such that the difference in travel time is sufficient for the pulses to 
arrive at distinct intervals in time. After the propagation losses and phase shifts are properly 
corrected, the incident and reflected pulses are deduced at the point of reflection. The ratio of the 
Fourier transforms of these pulses yields the reflection' coefficient and, with additional calculation, 
the impedance. Plane wave reflection is assumed at the surface since kR » 1. In cases of grazing 
incidence, the direct and reflected pulses may not have distinct arrivals due to the small difference 
in path length. Then, a second microphone is employed to measure the direct pulse at a point 
where the arrivals are still distinct. This measurement is extrapolated to the position of the 
original microphone. In the past, Davies and Mulholland (1979), and Cramond and Don (1984), 
have used methods similar to this one. 

A second method employed is referred to as the “Standing Wave Ratio” (SWR) method. 
The ratio of the Fourier transform of the direct plus reflected pulses to the Fourier transform of 
the extrapolated direct pulse yields an interference pattern oscillating about unity. The values 
at the maxima and minima of the pattern are used to approximate 1) the magnitude of the 
reflection coefficient at the mid-point between each maximum /mini mum pair, and 2) the phase 
of the reflection coefficient at each maximum and minimum. 

A third method employed is based upon an empirical model for impedance proposed by De- 
lany and Bazley (1970). This model relates the real and imaginary parts of the surface impedance 
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3- Typical results for determination of flow resistivity <r of carpet material on 
plywood. Shown are the predicted and measured SPL versus frequency curves 
for sound propagating over the flat plywood table top covered with a thin 
commercial carpet. In this context, the SPL refers to the sound level at the 
reflection point relative to the free field sound level extrapolated to that point. 
The source and receiver heights are 2 cm and 0.3 cm, respectively; the horizontal 
separation distance between source and receiver is 80 cm. From the data shown 
it is inferred that a is approximately 1600 rayls cgs. 
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to the ratio of frequency to an effective flow resistance parameter. This relation is expressed in 
the power law formula 

Re{Z : } = R^/PiC! = 1 + 9.08(//<r)-° 75 (5a) 

In»{Zj} = A'; j p\ Cj = ll-9(//<r)-°' 73 (56) 

where / is the frequency, p x c x is the characteristic impedance of air, is the specific acoustic 
impedance of the reflecting surface, and cr is an effective flow resistance parameter. The method 
consists of comparing analytical curves, based upon the theory of reflection of a spherical wave 
from an impedance plane, with experimental data of the SPL above the flat table top for particular 
source and receiver positions. The analytical curves are computed using an algorithm presented 
by Chien and Soroka (1975). The SPL is considered to be a function of the source and receiver 
positions, the characteristic impedance of air, the impedance of the surface, and frequency. Thus, 
the effective flow resistance parameter, cr, can be inferred from a comparison of analytical to 
experimental data once the other parameters are measured, calculated or picked. An initial 
result for carpet-on-plywood is shown in Fig. (3). Notice that despite the thinness of the carpet, 
ss 6 mm, the experimental curves are well fitted by analytical curves for certain values of cr. 
Similar results were reported previously by Embleton, Piercy, and Daigle (1983) as well as by 
Berry and Daigle (1987), all of whom used the same approach. 

Insertion Loss Measurements 

A set of experiments have been designed to measure the effects of the presence of the ridge 
on the sound field. These effects are expressed in terms of the insertion loss. In the following 
experiments, the insertion loss is taken as the ratio (in dB) of the Fourier transforms of the 
pressure at a point in the field with the ridge present to the pressure expected at that point with 
the ridge absent. The expected pressure at a point in the field is extrapolated from a reference 
pressure measured well above and in front of the ridge. 

The experiments are of a 2-dimensional nature in that all measurements are made in a vertical 
plane which is perpendicular to the axis of revolution of the cylindrical ridge and contains the 
point source. All of the experiments are performed for both the bare plywood and carpet-on- 
plywood surfaces and perhaps one other acoustically soft surface. The acoustic impedance of the 
surface of the ridge is assumed equivalent to that measured on the surface of the table for the 
same surface material or covering. 

Where possible, the incremental distance between successive measurements w r ill be on the 
order of the smallest length scale in the direction of measurement. The length scales 


L x = R/(kR) x f z 


L y = RI(kR)' lz 



NAG-1-566 Semiannual Report No. 6, page 77 


Proposal, page 17 

were deduced from the theory. At 40 kHz, L x and L y are approximately 25 cm and 2.5 cm, 
respectively. 

The first measurements were made of the insertion loss along the bare plywood surface of 
the ridge. The field microphone measured the combined direct and reflected pressure at 10 cm 
intervals over the span of the ridge. Additional measurements were made over two different 
vertical paths originating behind the ridge. The horizontal position of the first of these two paths 
was roughly 64 cm downstream from the apex of the ridge, while the second path was located 
roughly 210 cm downstream from the apex (Fig. 4). The measurements over each of the vertical 
paths were made at intervals of 1-2 cm. 

A similar set of experiments will be performed over the ridge with a carpeted surface, and 
perhaps one other soft surface, depending upon the results obtained with the carpet. In general, 
the measurements over each of the vertical paths will be made at intervals of 1-2 cm and will 
extend to heights at least twice that of the ridge or roughly 65 cm above the top surface of the 
table. This will be necessary in order to observe oscillations in the geometrical acoustic region at 
low frequencies 5 kHz). 

After all of the above measurements of insertion loss due to a single ridge have been con- 
ducted, a second ridge will be placed in series with the first and, experiments will be performed 
in order to deduce the impact of echoes from the second ridge upon the insertion loss as measured 
for the single ridge. The measurements themselves will be similar to those made for the single 
ridge. It is expected that echo effects near the back surface of the first ridge could be significant 
since in this region certain echoed arrivals from the second ridge may dominate rays diffracted 
from the first ridge. 


SUMMARY 

The objective of this proposed work is to produce experimental and numerical data portray- 
ing key aspects of high frequency acoustic diffraction over a curved surface. The experimental 
work consists of measuring the impedance of the curved surface as well as the insertion loss on, 
above, and behind the curved surface. Several different surfaces will be used in order to observe 
the effects of various surface impedances on the results. The data is to be presented in such a 
way as to describe well the transition from shadow to light. In addition, corresponding numerical 
data is to be produced from the theory as presented by Dr. A.D. Pierce. Hopefully, a comparison 
of the numerical data with the experimental data will yield insights into 1) the behavior of the 
diffracted field in the transition region and near a caustic such as that along the curved surface, 2) 
the importance of echoes relative to diffracted rays, and 3) the limitations of various theoretical 
assumptions such a hat of a locally reacting surface. 
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4. Plot of the insertion loss at approximately 209 cm behind the bare plywood 
ridge as a function of the dimensionless height parameter t/\ The point rj> = 0 
corresponds to a point on the shadow boundary (y = 0) for a given value 
of frequency and x. Positive (negative) values of V’ correspond to positive 
(negative) values of the height above the shadow boundary, y, for given values 
of frequency and x. The experimental data correspond to frequencies of 5 kHz 
(c), 15 kHz (o), and 25 kHz (•). The source was at a net distance of 243.7 cm 
from the apex of the ridge. The solid curve represents the insertion loss as 
predicted by the theory for knife-edge diffraction from a rigid surface. 
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ABSTRACTS OF TWO FORTHCOMING PRESENTATIONS 


The following pages reproduce two abstracts that have been submitted for presentation 
at two different conferences. The first one will be presented by Yves Berthelot at INTER- 
NOISE ’88, to be held in Avignon, France, August 30-September 1, 1988. The second one 
will be presented by James A. Kearns at the upcoming 115th meeting of the Acoustical 
Society of America, to be held in Seattle, WA, May 16-20, 1988. 
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INTERNOISE ’88: Y. H. Berthelot 

Les Effets de Diffraction Associds avec le Rayonnement d’une Source Sonore Proche 
d’une Surface ni Parfaitement Plane ni Parfaitement Rfeflechissante. Yves H. Berthelot, 
Allan D. Pierce, Ji-Xun Zhou, and James A. Kearns, School of Mechanical Engineering, Georgia 
Institute of Technology, Atlanta, Georgia, 30332-0405, USA. 


Lorsqu’une source sonore est proche d’une surface ni parfaitement plane ni parfaitement rigide, 
un champ sonore trfes complexe s’fetablit. Pour essayer de mieux comprendre cette classe de 
probifemes, des experiences de modfelisation sont en cours d’etudes pour verifier la theorie des 
expansions asymptotiques (Matched Asymptotic Expansions) appliqufees a la thfeorie de Fock-van 
der Pol-Bremmer pour la diffraction des ondes par une surface courbe d’impedance finie. En 
particulier, on prfesente des rfesultats (thfeoriques, numferiques, et experimentaux) modelisants le 
champ acoustique diffracte par une ou deux collines recouvertes soit par de la terre sfeche soit par 
de 1’herbe. On s’intferesse particuliferement k la rfegion de pfenombre acoustique et au phenomene 
d ’echos. [Recherche financfee par NASA-Langley]. 
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Yves Berthelot, Assistant Professor 
School of Mechanical Engineering 
Georgia Institute of Technology 
Atlanta , Georgia, 30332-0405, USA 
Tel : (404) 894-7482 
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Computational studies of the diffraction integral occuring in the MAE theory of sound propaga- 
tion over hills and valleys. James A. Kearns, Ji-xun Zhou, Yves H. Berthelot, and Allan D. Pierce 
(School of Mechanical Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332) 

An important class of problems pertaining to outdoor sound propagation is that of the 
diffraction which occurs when the ground is neither perfectly flat nor perfectly rigid. Such 
problems are encountered in the study of long range propagation of sound over hills and valleys. 
It has been shown previously [J. Acoust . Soc. Am., Suppl. 1, 79, S30-31] that the theory of 
matched asymptotic expansions allows one to express the diffracted field in terms of a complex 
integral involving Airy functions of complex argument. In some limiting cases, the diffraction 
integral reduces to some computationally very efficient forms: an equation based on geometrical 
acoustics in the the illuminated region, a creeping wave series in the shadow zone, and a knife-edge 
Fresnel diffraction integral far behind the ridge. In the present paper, the transition between these 
different regimes is investigated numerically by computing the general integral, and particular 
attention is given to the matching with the creeping waves series solution in the penumbra region. 
Computational results are compared with data obtained in laboratory scaled experiments. [Work 
supported by NASA Langley Research Center.] 
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If possible, I would like to see this paper presented in the same session as any paper by Gilles 
Daigle on outdoor sound propagation. 
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